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_ A convergent approximation is proposed for a mean field density-density correlation function in 

a system with a two-phase interface. It is based on a fourth-order expansion of the Hamiltonian in 
, ^ terms of fluctuations around the equilibrium profile. The approach is illustrated by one and three 

^ ^. dimensional calculations for systems characterized by the Ginzburg-Landau functional. 
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■ Density correlations in the vicinity of equilibrium planar interfaces have been extensively studied by numerous 
^ I authors Several methods, including the density functional theory |]l|,D, the capillary wave model ^j^, and 

the eigenstate expansion of the fluctuations were used to calculate the correlation function. Most existing 

approaches have to confront a problem of divergence, caused by the vanishing energy cost of rigid shifts of the 
interface in an infinite system without external potential. We will focus in this work on the Hamiltonian second 

•-^ ' derivative Sp^,i,^splz2) eigenstate expansion method (sometimes also called "field theoretical method", In 

the framework of this method the divergence manifests itself in the term (t>o{zi)(f>Q{z2) / Xq where (t>oiz) and Aq are the 

■ "zero" eigenstate and eigenvalue; <poiz) = dpQ{z)/dz where po{z) is the equilibrium density profile. The eigenvalue Aq 
M [ is zero for vanishing external localizing field. The traditional physical explanation for the divergence is the following: 

the density fluctuations corresponding to the interface shifts, to the first order of shift amplitude, have component 
only along the zero eigenstate, Ap{z) ~ po{z + Az) — pa{z) = dp{z)/dzAz + C'(Az)^, and free wandering of the 
interface as a whole results in the ambiguity in defining the density-density correlation function. In the framework of 
^ the capillary wave model ( ) this ambiguity is overcome by renormalizing the equilibrium density profile by taking 
into account the free wandering mode. As a result, average width of the renormalized interface diverges for vanishing 
external field, and the zero-order term contribution to the density-density correlation function tends to zero. However, 
after such renormalization, one loses information about instantenous and local fluctuations of the interface. 

In an attempt to resolve this ambiguity, we would like to focus our attention on the question of divergence of the 
zero eigenstate term The main motivation of our approach lies in the following. 

It is true that in an infinite system without a confining external field the energy cost of rigid shifts of the planar 
' two-phase interface is zero, Tl[p{zo + Az)] = Tl[p{zo)]. Yet this does not mean that the amplitude ao of the density 
fiuctuation proportional to the zero eigenstate, Sp{z) = aop'oiz) can grow infinitely without a free-energy penalty. 
Physically it follows from the fact, that although the position of the interface is undetermined, the actual values of the 
density near the interface cannot go significantly beyond the density of either of the bulk phases. For all the realistic 
Hamiltonian functionals describing stable systems, this non-divergence is controlled by a positive coefficient in front 
_ of the highest power of the density. For the Ginzburg-Landau functional 



^ . n[p] = J i{^^p{z)r + {1 - p'{z)r)dz (1) 

o'. ' 

O the latter is +p'^{z). The divergence that appears when only the harmonic (second-order) terms of the expansion 
^ , of the Hamiltonian are used is the direct consequence of neglect of all higher-order terms. Hence, to eliminate the 
; ' divergence of the zero-eigenstate term in the correlation function, it is natural to try using higher-order terms of 
rN [ the expansion of the Hamiltonian around the equilibrium density profile (see 0| where the equilibrium profile was 
• calculated using higher-order terms). It turns out that fourth-order terms are sufficient to keep the zero-eigenstate 
■ " " ' contribution finite. We propose an approximate method to consider these fourth order terms; this method is formally 
introduced in Section II. Concrete examples for the Ginzburg-Landau Hamiltonian for one- and three- dimensional 
systems are presented in Sections HI and IV. 

However, after eliminating the divergence of the zero-order eigenstate term, it is natural to ask, which eigenstate 
or combination of eigenstates of the second derivative matrix describe the macroscopic shifts of the interface Ap{z) — 
Pq{z + Az) — po{z)7 The expansion coefficient of this density fiuctuation along the zero eigenstate remains finite even 
for infinite shifts, 
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p+ca 

I Ap{z)p'Q{z)dz po(+oo)[po(+oo) - po(-oo)] 



(2) 



for Az +00. This statement could serve as another argument for a finite average value of the amplitude of 
fluctuation along Pq{z). The same is true for all other bound (localized) cigenstates. Yet the projection of the shift 
Pq{z + Az) — po{z) onto low- lying continuum states diverges as Az 00. For the Ginzburg-Landau Hamiltonian the 
mean- field equilibrium density profile is poiz) = tanh(z), and the first non- localized eigenstate of the second derivative 
matrix is proportional to (3tanh(z) — l)/2. The integral 

/+00 
(tanh(z + Az) - tanh(z))(3/2tanh(z) - l/2)dz 2Az (3) 
-00 

diverges linearly for Az — > 00. It means that since macroscopic shifts have no energy cost, the appropriate linear 
combination of low-lying continuum states with some of the coefficients growing proportionally to the magnitude of 
the shift also has no energy cost. The finiteness of the average values of all the continuum spectrum amplitudes 
is another artifact of second-order truncation; in particular, it is a result of the neglect of the mixing of different 
harmonics in the third and higher-order terms. Consequently, our approximation is convergent only because, after 
going to the higher-order expansion in the zero term, we stopped short of going beyond the harmonic approximation 
for terms containing the lowest-lying continuum states. 

However, the main contribution to the correlations near the interface comes from the bound states, while contin- 
uum states are more relevant for the correlations in the bulk phases. A possible merit of our approximation is the 
improvement in the accuracy of density correlation calculations in the vicinity of the interstate. To compare our re- 
sults to experimental data, theoretical calculations should be supplemented by a priori knowledge of the macroscopic 
localization of the interface. 

For a three-dimensional system the situation is similar. For square-gradient energy functionals, the perturbation 
Ap{x, y, z) = /Oo(z -|- f{x, y)) — po{z) of an initially flat interface has the energy cost 

[p'oiz)fdz / / \Vfix,y)\'dxdy (4) 

-C30 J —00 J —00 

For long- wavelength fluctuations, /(a:, y) ~ exp[i(fca;X-|- fcyj/)], |fc| ^ 1, the energy cost of such fluctuations vanishes as 
|fcp. However, considering the expansion of the density fluctuations over the system of eigenstates (now in 3D) of the 
Hamiltonian second derivative matrix, ^l)i{r) = l/L (f)i{z) exp{ikxx) exp{ikyy), the divergent contribution comes not 
from the terms with localized (jfi, but from the bottom of the continuum of z-coordinate eigenstates. To illustrate that, 
in section IV we calculate the convergent contribution to the correlation function from the (I)q{z) exp^ik^x) exp(ikyy) 
eigenstate. 



II. FORMALISM 



A density-density correlation function g(zi, Z2) is defined as a thermal average of a product of density fluctuations 
Ap{z) ~ p{z) — po{z) around the equilibrium density profile po{z): 

5(zi,Z2) = (Ap(zi)Ap(z2)). (5) 

For simplicity, in this section wc consider a one-dimensional case, p = p{z). Following we express it as a functional 
integral over all possible density profiles, 

(7(zi, Z2) - [l/Z] j Vp Ap(zi)Ap(z2) exp{~n[p]} (6) 

where Ti-lp] = H[p\/khT is a reduced Hamiltonian functional, and the partition function Z serves as the normalization 
constant, 

jvp exp{-n[p]} (7) 

In the mean-field approximation that we will use through this work, the equilibrium density profile pn{z) is deter- 
mined as the one that minimizes the Hamiltonian, 
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5p{z) 



= 0. 



(8) 



To evaluate (|6|) , we proceed by calculating the eigenstates of the integral operator with the kernel ^^^^^^^^^^^-^ 



P=PO 



4>i{z2)dz2 = \i4>i{zi). 

P=PO 



5p{zi)5p{z2) 

There is always a special eigenstate corresponding to 4>q{z) — dpo{z)/dz which has zero eigenvalue Aq = since 



(9) 



Sp{zi)Sp{z2) 



dpo{z2) , d SHlp] 
■ az2 = 



dZ2 



dz 5p{z) 



= 0. 

P=PO 



(10) 



The Hamiltonian is usually real and contains only even powers of differential operators; it makes the integral 
operator in Eq. Hermitian. The system of eigenstates (pi is complete and orthogonal; we also assume that it 
is normalized with unit weight function and all eigenfunctions (pi are made real. An arbitrary density fluctuation 
Ap(2;) = p{z) — po{z) can be expanded over the complete set of functions {(f)i{z)}: 



Ap{z)4>i{z)dz 



i=0 

Using the expansion the functional integral can be expressed as 

-22) = / ••■ / aiaj(f)i{zi)(l)j{z2) exp{-H[po + '^ak<t'k]} Y]_ 



(11) 



i=0 j=0 

where the normalization constant 



fc=0 



m— 



00 00 

. . . / exp{-'H[po + ^flfc^fc]} Y\_ 

•' k=0 m=0 



(12) 



(13) 



We assume here that the integrals in both numerator and denominator are convergent. As we mentioned in the 
Introduction, it is not true for the single specific direction in the space of coefficients {a^}, corresponding to the rigid 
macroscopic shifts of the interface. However, as a result of the approximations made below, this divergence will not 
affect the further calculations. 

Traditionally, Ti.[po + J^kLo '^kpk] is expanded to the second order around the equilibrium density profile, the 
orthogonality conditions for the are used, and the corresponding Gaussian integrals factorizcd As a result, 

the familiar expression for the density-density correlation function is recovered: 



{Z1,Z2) = ^{af)(j).i{zi)(l)^{z2) 



(14) 



i=0 



where (a|) = l/A^. As we already mentioned, the zero term diverges since Aq — 0. 

However, as we discussed in the Introduction, from physical considerations (ag) must have finite value. It is indeed 
the case if in Eqs. (|l2|), (|l3|) one goes to higher than second order in the expansion of 

In our case, the expansion of Ti.[p] up to the fourth order in the density fluctuation around the equilibrium profile 
is sufficient: 



Ti-ipa + ^ak(pk] « 

k=0 

II«M jM\ 't'k{z)dz + 



fe=0 



l5p(zi)(5p(22) 



_ ^k{z\)<pi{z2)Wdzj 

J = l 



P=PO 



(15) 
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if ^ akCLiam 

k.l,m—0 



s'n[p] 



Sp{zi)Sp{z2)Sp(z3) 

s^n[p] 



i'k{zi)(l>liz2)4>Tniz3)Y\_dZj 



Sp(zi)Sp{z2)Sp{z3)Sp{zi) 



P=Po 



The first order term is identically zero, the second-order terms are used in the traditional formalism, and the third 
and fourth-order terms are essential for our treatment. If we substitute this expansion into the expressions for the 
correlations function ([l^), (^3|), for non-pathological forms of the Hamiltonian functional 7i[p], it will produce a finite 
value for (ag). However, in their complete form Eqs. (|lj), (|l^), ( |l5| ) are hardly tractable. Assuming that we are far 
enough from a critical point, we use the following approximation to evaluate {a?^ We assume that the second-order 
expansion works well enough for all (of) — 1/Xi with i ^ and drop from Eq. (Uq) all terms that do not contain uq. 
In the remaining terms we replace all combinations of at, Uiaj, aiajak, {i,j,k} ^ by their average values obtained 
with the second-order expansion: 



(a^) = 0, 
(ciiaj) = 5ijj^, 

(aiUjak) = 0. 



(16) 



This approximations allow us to express (ag) in the form of the following integral: 

/ 2\ I exp{— ax"' — — jx — Sx'^}dx 



with coefficients a, (3, 7 given by 



00 



2^ A 

i=l 



J exp{— ax^ — Px^ — jx — Sx^}dx ' 



s^n[p] 



Sp{zi)Sp{z2)Sp{z3)Sp{z4,) 



^0 (^1 ) '/>0 (Z2 ) 00 (2:3 ) </>0 (Z4 ) d^^j , 



6p{zi)5p{z2)Sp{z3)dp{z4 

f f f s'n[p] 



t{zi)(j>i{z2)(t)o{z3)<j>o{z4) Y[d.Zj, 



5p{zi)Sp{z2)Sp{z3) 

s'n[p] 



{zi)(pt{z2)4io{z3) Y\_dzj, 



P=Po 



Sp{zi)5p{z2)Sp{z3) 



P=PO 



'o{zi)(j)o{z2)(poi,Z3) Yldzj. 

3=1 



(17) 



(18) 



III. GINZBURG-LANDAU HAMILTONIAN 



To illustrate our approach let us a consider a simple ID system with the Ginzburg-Landau Hamiltonian (|l|). The 
equilibrium density profile for symmetric boundary conditions p{—l) = —1, p(+0 = 1, ^ ^ 00 is pa{z) = tanh(z). The 
eigenvalue equation (||) takes the differential form 



- 2 — (j)^{z) - 4:<P^{z) + 12tanh^(z)0j(z) = A,0.,(z), 



and the third and fourth-order functional derivatives of 7i[/9] are: 

5'n[p] 



5p{zi)5p{z2)5p{z3) 



24tanh(zi)(5(zi - Z2)5{z2 - Z3) 



(19) 



(20) 
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Sp{zi)Sp{z2)Sp{z3)Sp{z4) 



24(5(zi - Z2)S{z2 - Z3)d{z3 - Z4). (21) 



In fact, the expansion ( [l5| ) is now exact for Ti-lp] since all higher-order variational derivatives are identically equal to 
zero. To proceed further we need to calculate a, 7, and 6 as defined in Eq. (p^). From parity consideration, the 
coefficients 7 and S are zero. 

The calculation of a is straightforward: 

Q r+°° dz 18 

9 / az ^22^ 



16 7-00 cosh«(^) 35 



The coefficient -j^ appears from the normalization condition: 



16 

f +00 



(/)^(z)dz = 1. (23) 



To evaluate /3 we need to know the eigenstates (t)i{z) of Eq. ( p^ 

9 - f+00 ^2(^) 



''^ '^~^J-oo cosh^(z) ^ 

The sum in Eq. ( p4[ ) contains a contribution from one bound state, {(t>i{z) = •\/3/2 sinh(z) cosh~^(z), Ai = 6) and 
continuum states (?> 2) To evaluate the sum over the continuum, we use expressions obtained in p] (Eq. (3.15 - 
3.20)): 

9 cj^^z) dz ^ 9 dk 1 1 ^fc(z) 

2^i_^ A. cosh4(z) 4i_^ 2^(P+4)2(fc2 + l)i-oo cosh4(z) ^ ^ 

where the continuum eigenstates (jikiz) are 

(j)k{z) = exp (iA:z)[l + fc2 + 3i/ctanh(z) - 3tanh^(z)]. (26) 



The normalization of 4)k{z) is taken into account in the fc-dependent factor of the outer integral in (|25|). Both 
integrations in (^) are straightforward and finally we obtain for (3 

/3 w ^ + 0.022 w 0.193 (27) 

The contribution from the continuum to the value of /3 is « 13%. We substitute these values of a and j3 into Eq. (|l 
and, in the first order in j3/^/a^ obtain for (oq): 

/ 2n ^ /a^^exp{-az'^}(l-/3z^)(j2: _ 1 ^^(P ~ ^fe^ _ ,-. . /9n^ 
/exp{-azn(l-/3^^)rf^ ' ^/^ - ^T^{\) ^ ^ ^ 

It is interesting to note that if one ignores all cross-terms ( contribution from j3) in Eq. (|l7|), (oq) « 0.471. For the 
Ginzburg-Landau Hamiltonian, « 0.27 plays the role of a small parameter in the approximations (|l6|) - (p^, as 

well as in (^8|). Hence there is a certain degree of numerical justification of heuristic assumptions made in ( p6| ) - 
(|l8|). To obtain the final result for the density-density correlation function (|l^ ) we use the following method. When 
the external potential V{z) = cz, linear in the coordinate z measured form the interface location z = 0, is added to 
the Hamiltonian, the zero eigenvalue becomes proportional to the coefficient in front of this term: Aq oc c. In |^ an 
expression for the correlation function Gc{zi, 22) (Eq- (4.17) in ||^) is obtained in the presence of an external potential 
F(z), with V{z) ctanh(z) when c ^ 0. It is shown that the zero eigenvalue Ag/c — > 1 for c — + 0. Since for i > \ 
all Ai go to constant limits when the external potential is turned off, the "truncated" correlation function g (without 
the zero eigenstate term) can be expressed as: 

5Zi,Z2 => r = Imi — cGc(zi,Z2). (29) 

^ — ' A, c^o dc 

i=i ' 
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To recover the "non-truncated" correlation function g{zi,Z2) we add to g{zi,Z2) from ( p9| ) the correct contribution 
form the zero-order term, 



4 cosh (zi) cosh (22) 



(30) 



with (og) given by (28). Sketches of the density-density correlation function ( |30| ) are presented in Figs. 1,2. It is 
straightforward to demonstrate that far from the interface {zi.Z2 3> 1), the density-density correlation function ( |30| ) 
goes to the correct bulk phase limit, 



exp(-2|zi - Z2I) 



(31) 




Fig. 1. A contour plot of the density-density correlation function g{zi,Z2). 
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Fig. 2. A 3D sketch of the density-density correlation function g(zi,Z2). 



IV. 3D CALCULATION 

The next logical step is to generalize this approach to a more realistic 3D system. For simplicity, we consider the 
same Ginzburg-Landau Haniiltonian as in (|l|) 

n[p]= li\Vp\' + {l-p^f)r)dr (32) 

The equilibrium density profile for the symmetric boundary conditions p(z = —I) = —1, p{z = +1) = 1, / — > cx) is the 
same as in the ID case, po{z) = tanh(z). The three-dimensional analog of the eigenvalue equation ([l^) reads 

- 2A^,(f) - A^,{r) + 12tanh2(z)?A,(f) = X^M^^, (33) 

Since the "potential" term in ( |3^ depends only on z, the variables can be separated, ipiir) = 4>i{z)^i{x, y). We assume 
that in the xy plane the system is confined to a square box of size L with periodic boundary conditions. Then the 
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xy components of the eigenstates of (^) are the normahzed plane waves £,j:{x,y) = {\ / L) eyi'p{ikxx) e'xp{ikyy) with 
^{x,y} ~ 27rn{a;_yj./L, n^^y-^ = 0,±1,±2.... Expanding density variations over the complete set of functions ip^j^ir) 
and taking the thermal average, we obtain for the density-density correlation function (compare to dTT 



i.ki j,k2 



exp{-H[po + E n d<^nk,- (34) 



Similarly to (13), Z is a normalization constant: 



m,ki n,k2 



Here and below we use the shorthand notation: for sums, products, and subscript indexes symbol k means the set 
of xy eigenstate indexes {nx,ny}. After expanding (34) to the second order in the density variation, an expression 



analogous to ( |14| ) is recovered: 

g{n,f2) = EKV'/';fc(^^i)^.fe(^^2), (36) 



i,k 



with (a^^) = l/('^i + 2fc2). As in the (Q), a similar problem arises for the Ao = eigenstate: for L ^ oo the sum on 

k diverges at the lower limit. 

To avoid the divergence occurring in the Ao = term, we suggest the same recipe as in the one-dimensional case: 
to continue the expansion of the Hamiltonian to the fourth order. For simplicity, we neglect mixing of the eigenstates 
with different Ai in the fourth-order term. Similarly to the ID case, where the relative contribution of mixing is given 
by the ratio of coefficients P/y/a « 0.27 (^ 27, 2|), the inclusion of this mixing here will not affect convergence but 



will slightly change the numerical values of the coefficients. For all A; ^ 0, the second-order result (|36| ) is sufficient, 
hence we perform fourth-order expansion only for the subspace of eigenstates with Aq = 0. The contribution from the 
Aq = eigenstate to the density-density correlation function can be expressed as: 

50 (ri , f2 ) = ^ E E / • ■ • / «ofci «o4 Kk, i^^Kk, ) 

fci k2 

exp{-W[po + E «0fe3^0fe3]} n^"ofe4' (3^) 

ka ki 

Zo = j ... j exp{-7^[po + E°ofci^ofeJ} H'^^ofc^- 

ki k2 

Using the orthogonality conditions for ^^(z,?;) and Eq. ( [l8| ) we obtain 



Zo 

ki 



where a is defined by (p2|). Direct evaluation of the functional integrals in (^8|) is impossible; however, a simple 
approximation will allow us to obtain a physically reasonable expression for go['ri,f2). The main contribution to 
the integral in ( ^ ) comes from a^j: with small |fc|, so in the first approximation it is natural to introduce an upper 
cutoff C in sums on k. In particular, we replace the infinite limits in all the sums and products on Ux and Uy in 
( ^ ) by the finite cutoff, \nx^y \ < VCL/2. It corresponds to a system-size-independent cutoff for a wavevector k with 
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\kx,y\ < VCtt. We select C in such a way that the term — 2 '^qJ^^i quadratic in a^^, can be neglected in the 
exponent, which decouples the integration over da„r: from the summation on nx,ny. Besides neglecting the quadratic 



term, we remove the functions i^nj^ir) from (|38|), and call the remaining expression A. 



JV N N 

j ■ ■■ j "fci'^^Pi-;^!!! «fe2]'} X{dak, 

^ ~ ~ N N ' 

with = CL^. This expression can be evaluated by using A^-dimensional spherical coordinates, = 



(39) 



A = 



1 



N 



(40) 



Here = /T{N /2) is the area of the N-dimensional sphere. By inspection, one can identify A as the sum of 
the first N terms of (a^), A = X^zlii'^^)- Therefore 



A 1/1 



Now we return to Eq. ( |38| ) and replace one sum in the fourth-order term in the exponential by 



After this substitution, Eq. (p8|) becomes a product of Gaussian integrals and its evaluation becomes trivial: 



\^ f ~2^'^'P[ii.kxX + kyy)] 
5o(fi,f2) ~ (t)o{zi)cj)o{z2)2_^a^ — , 



(41) 



(42) 



(43) 



where the "improved" average values a| are given by 



J al exp[-a|(2/c2 + V^/2)]dag l 
/ exp[-a2.(2fc2 + ^/^/2)]da^ ~ 4fc2 + 



(44) 



Taking the limit L — > oo and replacing the summation X^n^ n„ ™ (El) integration (L/27r)2/ J dk^dky, we 

obtain 



1 



kdk 



(471)-^ Jo + VaC/4 (47r)^ 2 



(45) 



Here Jo and Kq are Bessel and Modified Hankel functions of zero order and r±_ = \J {x\ — x^)^ + (yi — y-i)^ ■ For large 
positive q one has KQ{q) = Y^7r/2qexp(— g)[l + 0{l/q)]. Consequently, we can identify (aC)^/'*/2 with the previously 
introduced upper cutoff for the wavevector fc, i.e., with ^/Cn. It allows us to express the constant C through the 
known parameters of the system, C — aj (47r)"'. Finally, we write for the zero-eigenstate contribution to the correlation 
function: 



1 



ffo(n,r2) ~ (l)Q{zx)(l)Q{z2)J—^KQ\—^J{xx - X2Y + {v\ - V2?\- 
The contribution from other eigenstates, g{f\,r2\ is obtained by straightforward integration, 

5(^1,^2) = E'^j(^i)'^j(^2) ^^^''^^2 -^"[\/^^(^i ~ ^2)^ + {y\-V2Y\ 



(46) 



(47) 
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V. CONCLUSION 



We show that taking into account the fourth-order terms in the expansion of the Hamiltonian in the calculation of 
the density-density correlation function makes the previously divergent (for -D < 3) zcro-eigcnstatc term convergent. 
We also note that the macroscopic shifts of the interfacial profile are described not by the zero and other bound 
eigenstates, but by the combination of low-lying continuum eigenstates of the Hamiltonian second-derivative matrix. 
The inclusion of the convergent zero-order term allows us to improve the accuracy of the calculation of the correlation 
function in the vicinity of the interface. Our approach could be relevant for the experimental results obtained in the 
microgravity conditions using a local analytical method, e.g., scattering with a narrow beam, focused on the interface. 
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